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Abstract. I discuss the evolution of the redshift-space bispectrum via perturba¬ 
tion theory (PT) and large high-resolution numerical simulations. At large scales, 
we give the multipole expansion of the bispectrum in PT, which provides a natu¬ 
ral way to break the degeneracy between bias and Q present in measurements of the 
power spectrum distortions. At intermediate scales, we propose a simple phenomeno¬ 
logical model to take into account non-linear effects. N-body results show that at 
small scales the perturbative shape of the bispectrum monopole in redshift-space is 
preserved, breaking the hierarchical form valid in the absence of distortions. 


1 Introduction and PT Results 


The bispectrum, the three-point function of density fluctuations in Fourier 
space, is the lowest order statistic that carries information about the spatial 
coherence of large-scale structures. Non-linear perturbation theory predicts a 
characteristic dependence of the bispectrum on the shape of the triangle, which 
provides a signature of gravitational instability that can be used to probe the 
gaussianity of initial conditions and the bias of the galaxy distribution [^j. 
However, in order to use it in redshift surveys, redshift distortions due to 
peculiar velocities must be taken into account. In this talk, I discuss the 
evolution of the redshift-space bispectrum in PT and N-body simulations [Q . 

In the plane-parallel approximation, where the observed position s of a 
galaxy is given by s = x — f u z (x)z , with z a fixed direction, f(fl) « fP 6 , 
and u(a;) = —v(x)/(Ttf), with v(x) the peculiar velocity field; the density 
contrast in redshift-space reads [j7j 
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This equation describes the full non-linear density field in redshift-space in 
the plane-parallel approximation, and is the starting point for the perturba¬ 
tive approach in redshift-space. The term in square brackets describes the 
“squashing effect”, i.e. the increase in clustering amplitudes due to infall, and 
gives the standard Kaiser formula in linear PT ||. The exponential factor 
encodes the “fingers of god” (FOG) effect, which erases power due to velocity 
dispersion along the line of sight. From Eq. (jlj), it is straightforward to obtain 
the density field to any order in PT; in particular, the hierarchical amplitude 
Q s follows from second-order PT J|] 
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where the B s denotes the bispectrum and P°(fc) = b 2 ( 1 + |/3 + |/ 3 2 )P{k ) is 
the power spectrum monopole of galaxies in linear PT assuming deterministic 
biasing, and j3 = f /b. Decomposing into multipoles with respect to fi, where 
fj, = k\ ■ z, we obtain the tree-level monopole of the equilateral hierarchical 
amplitude [[fj: 

Qi% 0) = 5 A( 7 ,/3,6) /[986 (15 -|-10/3 + 3/3 2 ) 2 ], (3) 

where A = 2520 + 4410 7 + 3360/3 + 2940 7 /3 + 1260/3 2 + 441 7 /3 2 + 72/3 3 - 
636/3 3 — 146/3 4 , 7 = 62 /and b 2 is the non-linear bias. For no biasing this 
yields Q^l q = 0.464 for fi = 1. The quadrupole to monopole ratio of B eq is 

Rb = [5/(22A)](7392/3+6468 7 /3+3960/3 2 +1386 7 /3 2 + 264/3 3 -2316/3 3 -56&/3 4 ), 

(4) 

which for no biasing gives Rb = 11329/31394 = 0.36 for / = 1. Note that 
Eqs. (j^-Q) give a constraint on ( 6 , /?, 7 ) that is independent of the power spec¬ 
trum. These, together with Rp, the power spectrum quadrupole to monopole 
ratio statistic I 0 - can be inverted to obtain Cl, b and b 2 at large scales, 
independent of the initial conditions. 


2 Non-Linear Redshift Distortions and N-Body Results 

In order to describe the non-linear behavior of the redshift-space bispectrum, 
we introduce a phenomenological model to take into account the effects of 
velocity dispersion . For the power spectrum, we take p] 

P s (k) = b 2 P(k) (1 + /3m 2 ) 2 /[l + {k^a v f/2] 2 , (5) 

where a v is a free parameter that characterizes the velocity dispersion along 
the line of sight. It is straightforward to obtain the multipole moments of 
P s (k) in this simple model. Here we use the quadrupole to monopole ratio 
statistic, Rp, to fit a v , and then propose a similar ansatz to take into account 
non-linear distortions of the bispectrum: 

B s (ki,k 2 ,k 3 ) = B^ T (k 1 ,k 2 ,k 3 ) /[l +a 2 [(fci/ti ) 2 + {k 2 ^ 2 ) 2 + (fc 3 /z 3 ) 2 ] 2 <j 2 /2] 2 , 

( 6 ) 

where Bf T (ki,k 2 ,k 3 ) is the tree-level redshift-space bispectrum. Note that 
we introduced a constant a which reflects the configuration dependence of the 
triplet velocity dispersion. At z = 0, we take a ~ 2 for equilateral configura¬ 
tions, and a ~ 3 for k\/k 2 = 2 configurations, independent of cosmology. 



Figure 1 shows the results from this modeling and how it compares with 
N-body simulations, corresponding to 256 3 particles in a 240ft -1 Mpc box run 
by the Virgo Consortium. For the SCDM model, from Rp we obtain a v = 6 
(in units of Ho=100 ft km/s/Mpc) for z = 0 and a v = 2 for z=l (with a = 1). 
Using this we can predict Rp by taking multipoles in Eq. (||), the resulting Rp 
is shown as solid lines in Fig. 2, which despicts excellent agreement with the 
N-body simulation results. Similarly, for the ACDM model we obtain from 
fitting Rp, <j v = 5.5 for 2 = 0 and a v = 4 for z = 1 (with a = 1.75). The 
results for Rp are also in very good agreement with numerical simulations. 

The bottom panels in Fig. 1 show a comparison of simulations to the pre¬ 
dictions of PT and the model of Eq. (]g|) in fc 2 /fci = 2 configurations, for 
three different scales. We see in these panels that even though the tree-level 
PT prediction in real space (dotted) works reasonably well, the redshift-space 
counterpart (dashed) does not. On the other hand, the model in Eq. (Q) (solid) 
describes the N-body results very well. Even at the largest scale we probed, 
corresponding to a wavelength A ~ 50 Mpc/ft, the tree-level PT prediction in 
redshift-space would predict an effective bias b = 1.4, a quite significant dis¬ 
crepancy. This situation is similar to what happens with the power spectrum, 
redshift-space statistics are more affected by non-linearities than their real- 
space counterparts. In this respect, it is interesting to note that even in the 
linear dynamics, the exponential factor in Eq. (Q) can lead to a FOG effect at 
large scales that eventually makes Rp and I?b to become negative |J]. There¬ 
fore, the long range of the FOG effect seen in numerical simulations, should not 
be exclusively attributed to virialized clusters. These results strongly suggest 
the possibility of extending the leading-order PT results for the power spec¬ 
trum and bispectrum to smaller scales by treating the redshift-space mapping 
in Eq. (|Ij) exactly and approximating the dynamics using PT j?j]. The bottom 
right panel in Fig. 1 nicely illustrates the effect of cluster velocity dispersion on 
redshift-space correlations in the non-linear regime, whereas Q is very close to 
hierarchical, Q s has a strong configuration dependence. Eq. (|(j|) (solid) does an 
excellent job in predicting Q s , even at this considerable stage of non-linearity. 

Acknowledgements. This material is based on a paper in collaboration with H. 
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Supercomputing Consortium ( http://star-www.dur.ac.uk~frazerp/virgo/virgo.html ) 
using computers based at the Max Plank institut fur Astrophysik, Garching and the 
Edinburgh Parallel Computing Centre. 
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Figure 1: The top two panels show the quadrupole to monopole ratios i?p and 
i?B hr SCDM (top) and ACDM simulations. The full lines show the predictions 
of PT convolved with an exponential velocity dispersion model, Eqs. (|[](]). The 
bottom four panels show the power spectrum and the hierarchical amplitude 
Q in real (squares) and redshift (triangles) space for k 2 /ki configurations as a 
function of the angle 0 between k\ and k 2 for three different scales. 



















